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Abstract 


A combination of molecular quantum electrodynamics, perturbation theory, and ab 
initio calculations was used to create a computational methodology capable of esti¬ 
mating the rate of 3-body singlet upconversion in organic molecular assemblies. The 
approach was applied to quantify the conditions under which such relaxation rates, 
known as energy pooling, become meaningful for two test systems: stilbene-fluorescein 
and hexabenzocoronene-oligothiophene. Both exhibit low intra-molecular conversion 
but inter-molecular configurations exist in which pooling efficiency is at least 90% 
when placed in competition with more conventional relaxation pathways. For stilbene- 
fluorescein, the results are consistent with data generated in an earlier experimental 
investigation. Exercising these model systems facilitated the development of a set of 
design rules for the optimization of energy pooling. 
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INTRODUCTION 

The frequency of available light is often lower than desired for a given task. For instance, low- 
energy light can take advantage of an optical transparency window of biological tissue 7 and 
penetrate to ingested theranostic nanoparticles, 7 but the photoactivated release of cancer 
drugs requires energy in the UV range. 7 Within the same backdrop of internal medicine, 
the fluorescence of biomarkers is difficult to detect if it is in the same frequency as the light 
used to excited them. 7 The efficiency of photochemical processes in biological systems also 
suffers as a result of a mismatch of existing and desired frequencies of light. For instance, 
microalgae used in energy harvesting strategies do not efficiently utilize the red end of the 
solar spectrum for biomass production. 7 An analogous issue confronts those involved in the 
design of solar panels, since a large swath of the solar spectrum reaches earth at frequencies 
too low to be utilized by existing photovoltaic technology. 7 

Motivated by such disconnects between the desired light and that available, a number of 
strategies have been developed to upconvert light to higher frequencies. One of the simplest 
upconversion methods is the two-photon absorption (TPA) shown at left in Fig. [lj but 
this occurs only under extremely high intensity light. 7 7 Nanostructured condensed mat¬ 
ter, though, can be fabricated to independently absorb low-energy photons and subsequently 
combine the separate excitons into a single, higher energy state. Such upconversion processes 
are distinct from nonlinear optics in that the initial excitons can be created at distinctly sep¬ 
arate times and, typically, on differing nanostructures. These processes, for instance, allow 
the infrared spectrum to be utilized to carry out high-energy tasks in photovoltaics 7 7 7 , 
biofuel production and medical applications. 7 7 

In principle, upconversion can be carried out by sequentially exciting a single electron, 
known as Excited State Absorption (ESA), as shown at right in Fig. [lj However, it is 
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difficult to design systems in which ESA is efficient while ensuring that the resulting energy 
can be transferred before the exciton thermalizes. 7 Much more practical to implement is a 
two-electron process, shown in Fig. [2j in which a sensitizer molecule absorbs photons and 
transfers excitons to an activator. The latter are typically lanthanide atoms for which the 4f 
excitation is encapsulated within the 5s bonding orbital. 7 A low electron-phonon coupling 
results and the thermalization is consequently slow. 


Two-Photon Excited State 

Absorption (TPA) Absorption (ESA) 



(a) (b) 


Figure 1: Energy level diagrams for two-photon absorption (a) and excited state absorption 
(b). Times are indicated by t\ and £ 2 - Solid black lines indicate eigenstates. Yellow-black 
lines indicate virtual states. 

From the perspectives of biocompatibility, environmental safety, expense and resource 
abundance, it would be advantageous to use organic materials for such Energy Transfer 
Upconversions (ETU). High reorganization energy, and its manifestation as rapid phonon 
generation, has precluded any meaningful organic molecular ETU, though, because excited 
states thermalize into the first excited state faster than any other competing process. 7 It 
was this issue, in fact, that drove research towards lanthanide nanoparticles in the first place. 
This has motivated a consideration of three-electron processes 7 that completely side-step 
the thermalization issue. Energy Pooling (EP) is defined as a three electron process in which 
the energy of two electrons in excited states are transfered to the third electron. Such EP 
dynamics access virtual states to mediate upconversion, so there is no residence time during 
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which an excited state has the chance to lose energy to phonons. 

Energy Transfer Upconversion (ETU) 



Figure 2: Energy level diagrams for two-electron energy transfer upconversion (ETU) from 
donor molecules (D and D’) to an acceptor molecule (A). Solid black lines indicate eigen¬ 
states. Times are indicated by 


EP can occur via the two separate mechanisms illustrated in Fig. [3] The first is Accretive 
Energy Transfer (AET) in which one donor transfers its energy to a second donor, placing 
that donor into a virtual state. The second donor transfers its energy to the acceptor, and 
all time orderings contribute to the dynamics. The second pooling process is Cooperative 
Energy Transfer (CET) in which the energy of one donor is transferred to the acceptor plac¬ 
ing it into a virtual state. A second donor also transmits its energy to the acceptor, and the 
combined result is that the acceptor is in one of its eigenstates. As with AET, all time order¬ 
ing contribute the the total pooling rate. Although experimentally observed for lanthanide 


Accretive Energy Transfer (AET) 



(a) 


Cooperative Energy Transfer (CET) 



DAD' DAD' 


(b) 


Figure 3: Schematic depiction of two energy pooling (EP) mechanisms from donor molecules 
(D and D’) to an acceptor molecule (A). Solid black lines indicate eigenstates. Times are 
indicated by Yellow-black lines indicate virtual states. 
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nanoparticles, 7 7 the conditions under which EP will occur in organic condensed matter 
has not been explored. An earlier experimental study of a stilbene-fluorescein molecule 7 
suggests the possibility of organic upconversion providing an initial focus for the current 
computational investigation. Upconversion was observed to have occurred under high inten¬ 
sity excitation of the fluorescein donor resulting in excitation of the stilbene core. It was 
observed that unless two donors were bonded to each stilbene no upconversion occurred, 
eliminating TPA as a possible explanation. The paper concludes that either triplet-triplet 
annihilation (TTA) upconversion or some three-body resonance energy transfer (EP) could 
be responsible for the observed data but leans towards TTA. 

A quantum electrodynamics (QED) approach is combined with ab initio excited state 
structure calculations to quantify upconversion rates in assemblies of organic molecules. The 
codes are applied in a parametric study to formulate a set of design rules that can be used to 
guide the synthesis of systems that can efficiently carry out energy pooling. The computa¬ 
tional suite is subsequently applied to quantify upconversion processes in two condensed mat¬ 
ter assemblies: stilbene-fluorescein molecules as shown in Fig. |dj and hexabenzocoronene(HBC)- 
oligothiophene(OTP) molecules as shown in Fig. [5j 

The stilbene-fluorescein molecule studied is shown in Fig. |4} It is assumed that the 
molecule can be divided into a ^T-diaminostilbene core and fluorescein-5-isothiocyanate 
antennae components and each may be studied independently of each other. The interaction 
between the core and antennae are assumed to be dominated by dipole coupling. 



Figure 4: E-4,4 / -Di(fluorescein-5”-yl-thioureanyl)stilbene with two fluorescein donor arms 
connected to a stilbene acceptor. 
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Figure 5: Hexabenzocoronene (HBC) acceptor at symmetrically connected to six oligothio- 
phene (OTP) donor arms. 


METHODOLOGY 


Molecular QED provides a description of the competing energy transfer pathways, and those 
considered are shown in Fig. [6j Internal conversion and thermalization were not considered 
explicitly. Instead, it was assumed that electron-phonon coupling is sufficiently high that 
such processes, when available, dominate any relaxation dynamics. For instance, transitions 
to levels above the first excited state are assumed to immediately cool to the first excited 
state. 

All relaxation rates between initial and final states, |i) and |/), were derived from n tlci - 
order perturbation theory and so are expressed in terms of a coupling matrix element, W^ n \ 
and an energy density of states, 5: 



( 1 ) 


In turn, the coupling, W^ n \ has a functional dependence on the interaction Hamiltonian, 
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Figure 6: Schematic of energy transfer pathways considered in the computational paradigm: 
absorption (Abs), photoluminescence (PL), Energy Transfer Upconversion (ETU), Energy 
Pooling (EP), and Forster Resonance Energy Transfer (FRET). Internal conversion (IC), 
highlighted in red, is assumed to occur much faster than all other processes. 


V', that depends on the order of perturbation theory, n : 


W w 

= V ' 

( 2 ) 
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( 3 ) 
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( 4 ) 


It is often the case that the rates of two relaxation processes can be quantified using pertur¬ 
bation theory with each governed by the lowest non-zero interaction matrix element from 
Eqs. [2]-[4j In such situations, it might be tempting to conclude that the process of lowest 
order will be the fastest, but this is not always the case since molecular separation must also 
be taken into account. For instance, Forster Resonance Energy Transfer (FRET) can be 
more rapid than photoluminescence (PL) even though it is a second order perturbative pro¬ 
cess while PL is first order. Generally speaking, perturbative dynamics tend to have inverse 










polynomial relationships with separation distance, r, and the r-dependence of first, second, 
fourth, and sixth order processes are r°, r~ 6 , r~ 12 , and r~ 18 dependence, respectively. Fig. 
[7] shows that it is possible for each process to dominate relaxation over a range of separation 
distance, with closer distances associated with higher order processes in general. Scaling 
factors for the r°, r -6 , r -12 , and r -18 curves were chosen to be 1, 0.2, 2xl0 -3 , and 10 -7 , 
respectively. This choice was based on the assumption that the scaling factor should decrease 
exponentially with increasing order of perturbation theory. This suggests that energy pool¬ 
ing, a fourth-order perturbative process, will be observable only when molecular separations 
are very small. 



Figure 7: Idealization of perturbation processes of several orders as a function of molecular 
separation. The r°, r -6 , r~ 12 , and r -18 curves are scaled by factors of 1, 0.2, 2xl0~ 3 , and 
1CT 7 , respectively. 


For a given relaxation process, the interaction Hamiltonian can be expressed using the 
electric dipole approximation, and the form is the same for each of type of dynamics consid¬ 
ered in this work: 


E(r) 


V' = -qr ■ E(r) 


z £ ^^(k'AO^k'AOe^-) -^flt(k'A , )e ? (k'A')e-*( k, - r )], E k , 

k',A' 



( 5 ) 

( 6 ) 


In these equations, k is the wave vector, A is the polarization index, (kA) and a(kA) 
are photon creation and annihilation operators, £(kA) and ^(kA) are the unit polarization 
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vectors for each photon mode, c is the speed of light in vacuum, eo is the free space dielec¬ 
tric constant, and is the normalization volume. The rates of photoluminescence (PL), 
Forster Resonance Energy Transfer (FRET), and both pooling processes have previously 
been simplified into equations whose terms can be measured experimentally or calculated 
computationally: 7 777 


PpL = 


k 3^0a(D)y2 _ (E { - E { ) 3 {fJ? a W) 2 
3tT£qH 37T£o^ 4 C 3 


2tT | QfQi(D) T/ - /I -r*\ Bffii(A) 


rFRET = E 

k,l= 1 11 


S(E { -Ei) 


( 7 ) 
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with the electric-dipole coupling tensor, V and two-photon absorption tensor, a , given by 


D ^ftkR 


Vmn(k,H) = ^ 3 ((l±zkR)(S mn — 3R m Rn) - k R (Smn — RmRn)) (H) 


Mh 




pn pm 


(±fcl,±fc2) E ^ + r ,' TMl + B lc + !TT ficfe 


( 12 ) 


The following notation is used in these equations:' 


• k refers to the magnitude of the wave vector associated with a photon emitted from 
an excited donor; 

• /i xy represents the transition dipole moment from state y to state x on molecule z; 

• R, R/, and R" are the displacement vectors between molecules D and A, molecules D’ 
and A, and molecules D and D’ respectively. Hats indicate unit vectors while R , R ', 
and R" are vector magnitudes. 
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• 7 is a generic broadening taken to have a value of kpT (~25 meV). Varying this param¬ 
eter up to a maximum value of 100 meV caused less than an order of magnitude change 
in the rates. While this caused a quantitative change in the rates and efficiencies, there 
was no qualitative change in our results. This is consistent with other studies on the 
effect the broadening term on energy transfer processes between molecular systems. ? 

• C represents the possible virtual states that occur between the initial and final state. 


It is worth noting that an approximation is sometimes made to simplify the numerical 
evaluation of the two-photon absorption tensor, a, by reducing the sum over all states in Eq. 
?? to the sum of only the first and final state. However, this approximation is only reasonable 
under very specific conditions: (1) the molecule must posses a electric dipole moment that is 
larger than its transition dipole moments; (2) no singlet states must be closer in energy to the 
virtual state than the initial and final state; and (3) neither the initial or final state can have 
any excited states that are close in energy. The second condition makes this approximation 
invalid for the AET. For CET, this second condition can sometimes be met. However, for 
the acceptor molecules considered in this paper the electric dipole moment is essentially zero 
for CET so the first condition is not met. 

The square of the two-photon absorption tensor, otijOtM, is proportional to the rate of 
two-photon absorption, Tppq, given here in three different representations because each has 
a discipline in which it is commonly used: 


Ftp A 


1 2 0 

2 (hck) N/n 


I' 2 m _ 31 2 Im[ x ( 3 )] 

2 {hck) 2 * 7 An 2 eoHc 2 N/Q 


(13) 


Here P is the two-photon absorption coefficient, a( 2 ) is the two-photon absorption cross- 
section, and the imaginary part of the third-order electric susceptibility is ? 7 In 

addition, / is the irradiance of the light source and N is the number of molecules in the 
absorbing material. 

Summations are over the entire electronic eigenbasis although only a cutoff is imposed 
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in practice. The energy terms, E, refer to the electronic energy only and do not include 
the energy of photons that may be present in either the initial or final state. A Lorentzian 
function is used for the electronic density of states: 

1 7/2 

7 r (y/2) 2 + A E' 2 

The set of rate equations summarized above was populated using electronic structure 
analyses for the stilbene-fluorescein and hexabenzocoronene-oligothiophene geometries shown 
in Figs. |4] and [5j In the original experimental report on upconversion in stilbene-fluorescein, 
intermolecular structural information was absent except for a brief comment describing the 
material as a "red crystalline phase 1 '. 7 Because of this lack of structural information, the 
material was treated here as an assembly of interacting molecules without translational sym¬ 
metry. 

Molecular geometries were calculated using Density Functional Theory (DFT) and Time- 
Domain (TD)-DFT theory for ground states and excited states, respectively. The Def2- 
TZVP atomic basis set was used for all calculations. 7 However, a series of exchange- 
correlation functionals were considered that include: PBE, PBEO, B3LYP, and CAM-B3LYP. 7 7 7 7 
The impact on excited state properties of implementing the Tamm-Dancoff approximation 
(TDA) was also studied. 7 Absorption and emission spectra were produced directly from the 
TD-DFT calculation of the ground state geometry and excited state geometry respectively. 7 
For comparisons with experimental emission spectra, Kasha’s rule was enforced. 7 

Despite being referred to as ab initio calculations, an entire biom of density functionals 
have evolved for capturing electron exchange and correlation effects across the spectrum of 
condensed matter systems. Care must be taken in choosing a functional best suited to the 
system being studied, and this is particularly true for energy pooling. For instance, dis¬ 
crepancies in donor energy levels are doubly counted in the associated electronic density of 
states, and products of four transition dipole matrix elements make pooling rates especially 
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sensitive to errors in their estimation. As a guide for choosing exchange/correlation func¬ 
tionals, experimental absorption and emission spectra were measured and compared with 
those generated using PBEO, CAM-B3LYP, and B3LYP functionals. 7 777 

UV-vis absorption spectra for HBC was performed in solution using chloroform as a 
solvent. The emission data for oligothiophene was taken in a solution of chloroform. The 
oligothiophene samples were placed in 1 cm x 1 cm quartz cuvettes and excited with 380 
nm light. Emission data was taken over the range of 390 to 690 nm, measured 90 degrees 
relative to the excitation. These measurements were performed at room temperature using a 
PTI fluoronieter with an Ushio UXL-75XE xenon short arc lamp and a Hamamatsu R928P 
PMT tube operating at —1000 V, DC. 

As shown in Fig. [8j the CAM-B3LYP exchange/correlation functional was determined to 
be best for stilbene while PBEO gave the best match for fluorescein. The thermal broadenings 
used to generate the spectra were 250 meV for stilbene and 80 meV for fluorescein. In the 
case of stilbene, the experimental absorption data is for trans-stilbene rather than the trans- 
YY-diaminostilbene using for pooling calculations. A calculation was therefore performed for 
trans-stilbene in order to test the chemical accuracy of the exchange correlation functional 
with the functional chosen then used for all trans-^Y-diaminostilbene calculations. The 
results are shown in Fig. [8j The absorption and emission spectra for HBC and OTP, 
respectively, are shown in Fig. [9] The exchange/correlation functional that gave the best 
chemical accuracy were B3LYP and PBEO for HBC and OTP, respectively. The thermal 
broadenings used to generate the spectra were 100 meV for HBC and 150 meV for OTP. 
The secondary peaks shown in the experimental data of HBC are most likely the result of 
a combination of phonon and aggregation effects that are not accounted for in the model 
calculations. 
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Fluorescein Emission Spectra 



Stilbene Absorption Spectra 



Stilbene, experiment 

Stilbene, theory 
(TD-DFT, CAM-B3LYP) 

4,4'-diaminostilbene, 
theory (TD-DFT, 
CAM-B3LYP) 


300 350 400 

Wavelength (nm) 


450 


(a) Fluorescein-5-Isothiocyanate 


(b) Stilbene 


Figure 8: Comparison of spectra from TD-DFT calculations and experimental measurements 
using exchange/correlation functionals giving best match to data. These are identified in 
the figure legends. 

RESULTS 

Stilbene-Fluorescein 

Intermolecular and intramolecular dynamics were first considered for a single stilbene-fluorescein 
molecule, shown in its ground state configuration in Fig. |4^a). As a computational ex¬ 
pedient, the molecule was divided into a AT-diaminostilbene core and two fluorescein-5- 
isothiocyanate antennae with each moiety passivated and subjected to individual electronic 
structure analyses to obtain the required transition dipoles, electric-dipole coupling tensor 
and two-photon absorption tensor. The interactions between the core and antennae were 
assumed to be well-characterized by simple dipole coupling allowing Eqs. (7)-(10) to be used 
to calculate the rates of energy pooling and the competing processes. The rates are shown 
in Table [TJ The two orientations considered are displayed in Fig. 10 In simplest terms it is 
a comparison of intermolecular energy pooling and intramolecular energy pooling. 

In the intramolecular configuration, we see that the efficiency of the energy pooling is 
relatively small. The main factor for this is the center to center distances between the three 
molecules. Conversely, if we consider intermolecular system we see that smaller distance 
favors the energy pooling mechanism. There is an additional orientation factor that does 
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Oligothiophene Emission Spectra 



HBC Absorption Spectra 



(a) 


(b) 


Figure 9: Comparison of spectra from TD-DFT calculations and experimental measurements 
using exchange/correlation functionals giving best match to data. These are identified in 
the figure legends. 


Table 1: Rates and efficiencies of stilbene-fluorescein system. Where the sub¬ 


scripts a and b refer to the geometric orientations shown in Fig. 10 



r a (jus -i ) 

M%) 

r & (/rs i ) 

U%) 

PL 

8.1X10 1 

0.6 

8.1X10 1 

0.4 

GET 

1.9xlO u 

0.0 

3.2x 10 a 

16.9 

AET Total 

8.7xl0 -i 

0.0 

2.8x 10 1 

0.1 

ETU Total 

1.4xl0 4 

99.4 

1.6x 10 4 

82.7 


favor the intramolecular orientation, but it is small compared to the effect separation plays. 
The original experimental results did not quantify the efficiency of upconversion in such a 
way that a direct comparison could be made. ? 

The possibility of triplet-triplet annihilation (TTA) was also considered. In the geometry 
of the first excited state it was determined that the energy gap between the singlet state and 
the nearest lying triplet state was only 20 meV. This makes the possibility of intersystem 
crossing extremely likely. Additionally, when relaxed in the triplet state it was found that 
the energy gap widened to 400 meV, which means that the intersystem crossing in uni¬ 
directional; once the triplet state is formed it is very unlikely for it to revert back into a 
singlet. Since the rate intersystem crossing can not be calculated within our theoretical 
setting, the efficiency of this process is not known. 
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Fluorescein 


« 



Stilbene (Acceptor) 



Fluoresce 

(Donor) 



(a) Intramolecular 


(b) Intermolecular 


Figure 10: Molecular orientations of Stilbene-Fluorescein: (a) ground state geometry of a 
single molecule; (b) donor/acceptor orientations associated with optimized intermolecular 
pooling efficiency. S atoms are yellow, H atoms are white, C atoms are grey, O atoms are 
red, and N atoms are blue. The green (Stilbene) and pink (Fluorescein) shading is a guide 
to the eye and is also used to indicate the location of moieties not explicitly accounted for 
in the computational analyses. 

Hexabenzocoronene-Oligothiophene 

Hexabenzocoronene-oligothiophene (HBC-OTP) was studied next because the emission en¬ 
ergy of the OTP is over half the absorption energy of the HBC making energy pooling 
possible. In addition, the hexagonal shape of the HBC allowed for six antennae rather than 
two increasing the probability of two antennae being excited simultaneously, an important 
requirement for energy pooling. Last, it has been observed that HBC tends to form colli¬ 
mated liquid crystalline phases that can be aligned along a given direction by the application 
of an external electric held. This effect is believed to give rise to high carrier mobility. 7 
The geometry is shown in Fig. [5j Each component was assumed to be insulated and studied 
individually. The alkyl chains were methyl terminated for the HBC, but not for OTP due 
to the possible impact on the molecular geometry. 
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Fig. 11 shows four possible orientations for HBC/OTP system. This includes the simplest 
intramolecular case of the antennae being nearest neighbors, an intermolecular configuration 
of two OTP from neighboring molecules interweaving between the antennae of the HBC, a 
hybrid between the intramolecular and intermolecular configurations, and an intermolecular 
configuration in which the OTP are situation on either side of the HBC at a distance of 7A. 

Table [2] gives the results of the rates and efficiencies for each of the configurations of 
the HBC-OTP system. The efficiency of energy pooling in the intramolecular system is 
only about 0.4%; high enough to be experimentally measurable but considerably less than 
desired. The main cause of this low rate is the relatively large distance between the geometric 
centers of the core and antennae. The fourth configuration is by far the most efficient in 
any case. This is due to the small center-to-center distance of the core and antennae in 
this configuration. If the distance factor is changed to match the distance in the first three 
configurations, the fourth configuration would be less efficient than the others due to the 
direction of the transition dipole moments. 


Table 2: Rates and efficiencies of hexabenzocoronene-oligothiophene system. 
Where the subscripts a, b, c and d refer to the geometric orientations shown in 

Fig. [TT| 



r a (/rs i ) 

0a(%) 

r b (/is i ) 

M%) 

r c (/is i ) 

M%) 

r d(/ws i ) 

M%) 

PL 

6.8x ICE 

5.5 

6.8x 10 2 

0.5 

6.8x lO 2 

1.6 

6.8x 10 2 

0.0 

CET 

5.5x 10 1 

0.4 

8.5x 10 3 

6.6 

6.7xl0 2 

1.6 

3.5x 10 y 

99.0 

AET (D-D’) 

1.6xl0 -2 

0.0 

1.9xlO u 

0.0 

2.0x 10 u 

0.0 

1.3xl0 3 

0.0 

AET (D’-D) 

1.3xl0" 2 

0.0 

3.1xlO u 

0.0 

8.7x10“^ 

0.0 

1.3x 10 a 

0.0 

ETU (D-D’) 

5.9x lCF 

48.2 

5.9x 10 4 

45.9 

1.5xl0 4 

33.7 

1.8xhP 

0.5 

ETU (D’-D) 

5.6xlCP 

45.8 

6.1xl0 4 

47.0 

2.7xl0 4 

63.2 

1.8x 10 a 

0.5 


CONCLUSIONS 


Organic molecular assemblies tend to have strong electron-phonon coupling and a concomi¬ 
tant rapid internal conversion of excess exciton energy into heat. This is a primary obstacle 
in the design of organic upconversion materials. In such settings, energy pooling offers the 
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(a) Single HBC-OTP molecule 



OTP (Donor) 


(b) Co-planar, inter-molecular geometry 



(c) Mixed inter-molecular, intramolecular geome¬ 
try 



HBC 

(Acceptor) 


OTP (Donor) 

(d) Co-facial, intermolecular geometry 


Figure 11: Four representative geometries for HBC-OTP. S atoms are yellow, H atoms are 
white, C atoms are grey, and N atoms are blue. The green (HBC) and pink (OTP) shading 
is a guide to the eye and is also used to indicate the location of moieties not explicitly 
accounted for in the computational analyses. 
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prospect of efficiently up-converting excitonic energy by avoiding states in which hot exci- 
tons can lose energy to phonons. However, the conditions under which such relaxation rates 
become meaningful is not understood, motivating the current investigation. A combination 
of molecular quantum electrodynamics, perturbation theory, and electronic structure cal¬ 
culations was used to develop a simulator that calculates energy pooling rates along with 
other relaxation processes. These include photoluminescence and energy transfer upconver- 
sion but not triplet-triplet upconversion. Internal conversion was assumed to be much more 
rapid than any other relaxation dynamics. An important computational expedient was to 
assume that the interaction between donor and acceptor units is well-approximated by elec¬ 
tric dipole coupling. This allowed rate expressions to be developed that can be populated 
with electronic structure data associated with isolated units-e.g. transition dipole elements 
for isolated donor molecules and isolated acceptor molecules. 

Two specific example systems were chosen to exercise the methodology: stilbene-fluorescein 
and hexabenzocoronene-oligothiophene (HBC-OTP). In the case of the stilbene-fluorescein, 
it was found that intra-molecular energy pooling is measurable but with an efficiency less 
than 0.0%. This can dramatically improved, however, to approximately 17% efficiency for 
geometry-optimized inter -molecular dynamics. An examination of excited singlet and triplet 
energy levels, though, indicates that there is only small energy gap between the first singlet 
and the nearest triplet state. This leads us to conclude that triplet-triplet annihilation is 
also a reasonable interpretation for the upconversion observed in previous experiments, but 
the computational estimate of triplet-triplet annihilation is beyond the scope of the present 
investigation. For HBC-OTP, the efficiency of intramolecular energy pooling is estimated 
to be 0.4%, primarily limited by the large spacing between donor and acceptor components. 
On the other hand, an optimized co-facial orientation that could be associated with either 
inter-molecular dynamics or OTP arms bent back onto the HBC yields pooling efficiencies 
of 99%. Cooperative Energy Pooling was found to be much more efficient that Accretive 
Energy Pooling in all cases. We conclude that the future design efforts should focus on 
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inter-molecular dynamics that employ the following design guidance: 


• Energy pooling has an extremely sensitive r -12 distance dependence as compared with 
compared the r~ 6 dependence of energy transfer upconversion. This suggests that 
donor-acceptor distance should be minimized to promote energy pooling. As separation 
between components is reduced, though, multipole effects will become increasingly 
important as will Dexter processes. 

• Both donor and acceptor components should be chosen to have a large Two-Photon 
Absorption (TPA) cross section. This is because both pooling rates correlate with 
the TPA tensor which are directly related to the TPA cross-section (see Methodology 
section). As a result, the significant body of nonlinear optics data on TPA materials 
can be fruitfully applied to exciton dynamics as well. 

• Donor molecules should be used that cannot sequentially absorb two photons i.e. can¬ 
not carry out Energy Transfer upconversion (ETU). Along the same lines, energy 
transfer from acceptor to donor should be designed out. In both cases the intent is to 
exclude the possibility of hot excitons losing energy to heat. 

• All upconversion rates depend strongly on the relative orientation of donor and acceptor 
components. Our simulations show that it is possible to lower the coupling between 
the donors while maintaining a strong coupling between the donors and the acceptor. 
This would decrease the rate of ETU so as to reduce the subsequent generation of heat. 

ACKNOWLEDGMENTS 

This research is supported by the NSF SOLAR Grant CHE-1125937. All computations 
were carried out at the Golden Energy Computing Organization, Colorado School of Mines. 
We are pleased to acknowledge useful discussions with Prof. Gregory Scholes, Princeton 
University. 


20 



